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Energy Relaxation at a Hot-Electron Vortex Instability 
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At high dissipation levels, vortex motion in a superconducting film has been observed to become 
unstable at a certain critical vortex velocity v* . At substrate temperatures substantially below Tc, 
the observed behavior can be accounted for by a model in which the electrons reach an elevated 
temperature relative to the phonons and the substrate. Here we examine the underlying assumptions 
concerning energy flow and relaxation times in this model. A calculation of the rate of energy transfer 
from the electron gas to the lattice finds that at the instability, the electronic temperature reaches 
a very high value close to the critical temperature. Our calculated energy relaxation times are 
consistent with those deduced from the experiments. We also estimate the phonon mean free path 
and assess its effect on the flow of energy in the film. 

PACS numbers: 71.10.Ca, 71.38.-k, 72.10.Di, 72.15.Lh, 73.50.Fq, 74.25. Fy, 74.72.Bk, 74.78.Bz 



I. INTRODUCTION 

When a film of a type II superconductor is placed in 
a magnetic field large enough to permit penetration of 
vortices, a transport current in the film acts on the vor- 
tices through a Lorentz force that is opposed by a pin- 
ning force and, eventually, by a drag force. When the 
Lorentz force exceeds the pinning force, the vortices are 
set into motion and the drag force comes into play. When 
the Lorentz force is substantially larger than the pinning 
forces but the transport current is still small compared to 
the depairing current, previous experiments^*^ showed 
that the resulting dissipation is reasonably well described 
by the Bardeen-Stephcn (BS) modeli In this region it is 
Ohmic, but as the current is increased, it becomes non- 
linear and eventually reaches an instability manifested 
by a discontinuous increase in voltage. At temperatures 
not far below the critical temperature, the instability has 
been studied in a classic paper— by Larkin and Ovchin- 
nikov (LO). They showed that the electron distribution 
departs from a thermal distribution at high vortex veloc- 
ities, changing the superconducting order parameter and 
altering the drag force on the vortices. They predicted a 
non-linearity in the current- voltage characteristic and an 
instability in the vortex motion when the vortices reach 
a critical velocity v* . The LO instability is due to a de- 
crease in the drag force with increasing vortex velocity, 
accompanied by a decrease in vortex size. LO showed 
that the critical velocity is independent of the magnetic 
field. Early experiments on low-Tc systems^ confirmed 
Larkin and Ovchinnikov's results and predictions. Subse- 
quent experiments on YiBa2Cu307_5 (YBCO) by Doet- 
tinger, Huebener, Gerdemann, Kuhle, Anders, Trauble, 
and Villegier— and by Xiao and ZiemannS, also confirmed 
LO behavior. 

However, experiments carried out at lower 



temperature; 



,9,10 



on YBCO, showed a non-linearity 



and instability with a very different dependence of v* on 
the magnetic field B. Analysis 9 ^ showed that the new 
behavior could be accounted for by a simple model in 
which the electron gas has a thermal-like distribution 



function characterized by a higher temperature than 
the lattice and bath. Larkin and Ovchinnikov did, in 
fact, suggest this possibility in their original paper- 
without exploring its consequences. As the electron 
temperature rises, the resulting increase in resistivity 
causes a decrease in current above a certain electric 
field and hence a non-monotonic response. This model 
yields a critical vortex velocity v* at instability that is 
proportional to l/vB, as seen in the low-T experiments. 
Some of the essential consequences of such a hot-electron 
instability were calculated in our earlier papers and 
shown to be consistent with experimental observations. 

In the present work some of the simplifying assump- 
tions and restrictions in the previous calculations have 
been removed and more complete calculations have been 
carried out: 

1. ) The rate r" 1 of transfer of energy from the elec- 
tron gas to the lattice — which plays a crucial role in de- 
termining the electron temperature — was taken as a con- 
stant in previous discussions of the model. In this paper 
we show that it can be expected to have a strong tem- 
perature dependence. This temperature dependence of 
r e is now included in our numerical calculations of the 
current-voltage curves. We find that the general shape 
of the current-voltage relation is not very sensitive to 
the temperature variation of t £ because the electron gas 
passes rapidly from the bath temperature to a tempera- 
ture not far below Tc before any significant non-linearity 
is manifested. This is a consequence of the very small 
low-temperature specific heat of a superconducting elec- 
tron gas. However, the strong temperature variation of 
the relaxation time gives a sensitive measure of the elec- 
tron temperature. Evaluation of r e from the data near 
the instability point indicates an electron temperature 
much higher than the bath temperature, supporting the 
heated electron picture of the instability. The calculation 
of this electron-lattice energy relaxation time is presented 
in Section III below. 

2. ) In our previous work 9,10 , we assumed the film 
thickness to be negligible compared to the phonon mean 
free path, so that the phonon temperature is uniform 



2 



throughout the film. In this work, we remove this as- 
sumption and take phonon lifetime effects into account. 
The phonons will not necessarily be at the bath tem- 
perature, and will have a non-thermal distribution which 
varies with position in the film. Phonon lifetime effects 
can be taken into account following work by Bezuglij and 
Shklovskyi^, who solved the phonon kinetic equation for 
a thin film. The non-thermal phonon distribution found 
in this solution can be incorporated into our calculation 
of the energy transfer rate, and provides a correction to 
our earlier results. This result is derived in Section IV 
below. 

We begin in Section II by giving a description of the 
model presenting some new numerical results for the 
current-voltage curves under various conditions and for 
the critical parameters at the instability. 



II. MODEL FOR INSTABILITY 

The macroscopic fields in a type II superconductor car- 
rying a transport current are related to the velocity of the 
vortices by the fundamental relation 



E 

v = -c, 



(1) 



which follows from the law of induction. This equation 
can be used to find the electric field once the vortex ve- 
locity is determined by considering the fundamental dis- 
sipative processes in the medium. Elastic forces can be 
shown to be negligible. One of the dissipative processes is 
the scattering of normal electrons in the vortex core and 
quasiparticles outside the core first treated by Bardeen 
and Stephen. They found that the transport current den- 
sity j is expressed in terms of the upper critical field H C 2 
and the normal resistivity p n by 



. _ H C 2 E 
Pn B' 



(2) 



Later treatments^i^i^ take into account the relaxation 
of the order parameter during passage of the vortex first 
treated by Tinkham 12 . They give results which vary with 
circumstances, but agree with Eq.J2J) to within a numer- 
ical factor of order 1 . 

These energy dissipation mechanisms raise the energy 
of the electrons, and this energy subsequently relaxes to 
the lattice. The assumption of our model is that the 
electron-electron scattering time is small enough com- 
pared to the electron-phonon inelastic scattering time 
that the electron gas remains in internal thermal equi- 
librium at a temperature higher than the lattice temper- 
ature. The plausibility of the assumption can be checked 
by estimating the cross-over temperature below which 
electron-electron scattering is dominant. The standard 
estimates^ of the scattering rates = tjcf/T 2 and of 
Tgp 1 = rfuj'jj /T 3 then give a cross-over temperature of the 
order of 100 K for parameters appropriate to YBCO 25 . 



This temperature is indeed higher than the range of in- 
terest in the experiments. 

Changes in the energy density of the electron gas can 
be described by a rate equation that includes the work 
done by the electric field and the exchange of energy with 
the lattice. If we assume that the exchange can be de- 
scribed approximately by an energy relaxation time r e , 
then the equation is 



du _ _ u(T') - u{T p ) 
dt~ J T e (T',T p ) ' 



(3) 



where r e can depend on the phonon temperature T p as 
well as on the elevated electron temperature T". We ar- 
gue below that the dependence of r e on T p is weak enough 
to be ignored in the relevant range of temperatures and 
the relevant energy transfer rates between the lattice and 
the bath. The quasiparticles transfer the energy they re- 
ceive from the transport current to the lattice at a rate 
much higher than it is radiated back, and the energy 
then flows from the lattice to the bath. Thus r c can be 
assumed to depend only on T", and we can write the 
steady-state equation 



jEr e = I c(T)dT, 

IT* 



(4) 



where the energy difference in Eq. @ has been expressed 
in terms of the electronic specific heat per unit volume. 

Equations (J2J - (J3J determine the relationship between 
the electric field, the current density and the tempera- 
ture. The temperature dependence of the specific heat 
and the upper critical field are taken from standard BCS 
theoryii. In calculating the specific heat, the tempera- 
ture dependence of the gap was taken from BCS theory 
and its magnitude was multiplied by a factor to give the 
observed zero-temperature gap 18 and critical tempera- 
ture. In the next section, we calculate the energy relax- 
ation time and its temperature dependence. 

Typical results of the model are presented in the follow- 
ing figures. Fig. ^ shows the calculated current-voltage 
curves for different magnetic fields. The shape of the 
curves is in general qualitative agreement with the ex- 
perimental data shown in Fig|2 The strength parameter 
b of the electron-phonon coupling, defined in the next 
section, was adjusted to obtain agreement with the val- 
ues of E and j at the peak. We comment on the choice 
of b in the next section. 

The onset of the unstable region in the current-voltage 
response does not require explicitly invoking the forces 
on the vortices in treatment of the model. Rather, the 
instability appears in the result as a region of negative 
differential conductivity, where j decreases as a function 
of E. The region begins at the value E* of the field 
that can be determined by calculating dj jdE from Eq.Q, 
setting the result equal to zero, and solving for E: 



E* = 



C Pn B 



H c 2T e + {H c 2)'TeY 



(5) 
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where primes indicate differentiation with respect to 
temperature^i The experimentally well-verified \f~B de- 
pendence of E* follows provided the temperature T* at 
the instability is independent or weakly dependent on B 
so that the temperature-dependent factors C, r e , and 
H C 2 in Eq.J^J) remain independent of B. This result is 
a consequence of our model, since we have explicitly ex- 
cluded a field dependence for these quantities and taken 
p n to be temperature and field independent. Although 
VolovilsiS has shown that the specific heat has a in- 
dependence in type II materials above the lower critical 
field, we have checked that his scaling prediction at low 
temperatures gives only a weak dependence in the range 
of fields B « H C 2 relevant to our experiment. Fig. El 
shows the change in the electron temperature as a func- 
tion of the applied electric field. The rise in temperature 
and corresponding decrease in H C 2 result in decreasing 
differential conductivity which leads to the instability. 
Fig. 0] shows the relatively small effect of increasing the 
phonon temperature, up to about 40K, on the final tem- 
perature T* reached by the electron gas. 
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FIG. 2: Experimental curves of current density vs electric field 
in YBCO for flux density (beginning with the upper curve) 
B = 3, 5, 8, 11, 14, and 16 Tesla. 



III. ENERGY TRANSFER RATE 

The total rate at which energy is radiated by the 
heated quasiparticle gas to the lattice can be calculated 
by standard mcthodsSS^. The two contributing pro- 
cesses, phonon emission and quasiparticle recombination 
with emission of a phonon, are illustrated schematically 
in Fig. In the following calculation, applicable to d- 
wave superconductors, we assume the Fermi surface to 
be a cylinder of radius kp and height 2tt/cq, cq being the 
c-axis lattice constant. The rate for emission of a phonon 
of momentum q = k— k' by a quasiparticle of momentum 



k is 

W = J^Kfh J J ^^MM^Siu-uj^E^-E^-huj). 

(6) 

This rate can be expressed in terms of the electron- 
phonon spectral function a 2 (u>) F (ui) , defined by 

" 2 M^M= (2^3/? / t?! Mk - k 'l 2 ^- a;k - k ')' W 

where dS' is an area element on the Fermi surface and 
v' F is the Fermi velocity. If the quantity | A^k— k'l 2 ^^ — 
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FIG. 1: Effect of increasing flux density on the current density 
vs electric field curve calculated in the model with variable r e . 
Values of B beginning at the upper curve are 3, 5, 8, 11, and 
14 Tesla. 
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FIG. 3: Calculated electron temperature T' vs electric field 
E for B = 3,8, andl4 Tesla, with B increasing from left to 
right. 
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^k-k') is replaced by its average over the Fermi surface, 

2ircoh 2 VF 2 



I y ^'ix k _ k ,i 2 ^- Wk -kO 



the transition rate becomes 



Vk F 



-a 2 F(u), 



f d 2 k'a 2 F{EK-Ev)/K), 
3c h-y J 



(S 



(9 



where d 2 fc' = k'dk'dO' and 0' is the azimuthal angle o 
the Fermi cylinder. The total energy transfer rate wa 
obtained by integrating w over initial quasiparticle enei 
gies, and reduces to 



P, = 



u ^3j- J d 2 k J d 2 k'a 2 F(E-E')/h)(E~E')g e (l 

(1C 

The factor g e (E, E') contains the occupation factors an 
coherence factors 



g e (E,E')=f(E)(l-f(E>))[(l 



AA' cos 26» cos 26" 
EE' 

for the initial and final states. The second term in the 
coherence factor integrates to zero because of the d-wave 
symmetry of the order parameter. After transforming 



the momentum integrals into integrals over the quasipar- 
ticle energy and performing the azimuthal integrals, the 
energy transfer rate reduces to 




FIG. 4: Electron temperature at the instability vs. phonon 
temperature 
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J EdED(E, A) J E'dE'D(E',A)F((E - E')/H)(E - E')g e (E,E'). 



(12) 



Here D(E, A) is the d-wave density of states 



dB 



VE 2 - A 2 cos 2 20 



±K(P 2 

A V , 



±K(&), E> A, 



(13) 



where K denotes the elliptic function of the first kind. 
In this way of calculating, there is some averaging over 
the Fermi surface, but the characteristic d-wave density 
of states with its logarithmic singularity at E = A has 
been retained. 

The spectral function a 2 F(oj) is assumed to be of the 
form bu> 2 appropriate for acoustic phonons. Although op- 
tical phonons are present in high temperature supercon- 
ductors, they are assumed to make only a negligible con- 
tribution to the thermal conductivity responsible for car- 
rying energy from the heated electrons to the bath. The 
spectral function is cut off at the Debye frequency ujd- 
The constant b measures the strength of the electron- 
phonon interaction. This is usually expressed through 
the electron-phonon coupling constant A, defined as the 
integral of 2a 2 F(uj)/lu over frequency, which has a value 
of order unity in most superconductors. Since we are 
only considering acoustic phonons, this provides an up- 
per limit for b of about lCP 3 meV~ 2 . This is consistent 





FIG. 5: Diagrams contributing to the energy transfer rate 
from the quasiparticle gas to the lattice. 



with the magnitude of the electron-phonon matrix ele- 
ment quoted in Ref.|2l|. which gives b w 6 x 10~ 4 . It is 
also in the same range as the values extracted from neu- 
tron scattering data by Kaplan et alH2> for low temper- 
ature superconductors. The deformation potential ap- 
proximation as cited in Ref.[23, on the other hand, gives 
a value of the order of 3 x 10~ 7 mcV~ 2 . The choice that 
gives the best peak values for the peak current and for E* 
is 5.2 x 10~ 6 meV~ 2 , more in line with the latter value. 
The curves shown have taken this best fit value for b. 
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Introducing the dimensionless variables x = E/T, y — 
E'/T, and z — A/T, the energy transfer rate takes the 
form 



P e = GVT 5 $ f 



with 



G = 



3 7 6 



47T 3 fi 3 ' 

and the dimensionless integral $ e given by 

D(x,z)D(y,z)xy(x ~ yf 



d.r 



e x + e x-y + e -y + 1 



(14) 



(15) 



(16) 



D(x, y) is the dimensionless form of the density of states 
Eq. ^21- The quasiparticle recombination process gives 
an expression of the same form with <I> e replaced by 



d 

dx 
Jo 



d x d D(x,z)D(y,z)xy(x + y) 3 

y e x + p x+y + e y + I ■ V J 



The value of the specific heat constant 7 is taken from 
measurements of W. C. Lee, et. as 0.0331 meV -1 

nm~ 3 , giving G — 1.257 x 10 6 meV -1 nm~ 3 s" 1 when b 
is normalized to give A = 1. 
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FIG. 6: Calculated energy transfer rate in meV/nm -s of elec- 
tron gas at temperature T to a zero temperature lattice. The 
lower curve represents the phonon emission process and the 
upper curve represents the quasiparticle annihilation process. 



Fig. shows the calculated energy transfer rates for 
emission and recombination. Fig. [7| shows the energy 
relaxation time r c found by equating the total transfer 
rate from both processes to the last term in the rate 
equation (Eq. 0). 

We note the following properties of the energy transfer 
rate: 



1.) The rate of energy transfer from the electrons 
to the lattice at any given temperature is equal to the 
rate of transfer from the lattice to the electrons at the 
same temperature. Indeed, the rate from phonon emis- 
sion balances the rate from absorDtion and the rate from 




10 20 40 60 80 
T(K) 

FIG. 7: Energy relaxation time r E as a function of tempera- 
ture. 



quasiparticle recombination balances the rate from pair 
creation. These results can be demonstrated in the de- 
formation potential approximation, where the matrix el- 
ement for phonon emission and absorption depends only 
on the phonon energy v. For example, the rates for emis- 
sion and absorption and for quasiparticle recombination 
and creation can be written 



P e = G J (n(u) + l)f e (E,E - v)N{E)N(v ~ E)v z dEdv 
P a = G J n{v)f a {E,E + u)N{E)N{E + v)v 3 dEdv 1 
P r = G J \n{v) + l)f T (E,v-E)N(E)N{v- E)u 3 dEdv 
P c = G I n(u)f c (E,v - E)N(E)N{v - E)v 3 dEdv, 
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where n(y) is the phonon occupation number at the given 
temperature and N(E) is the quasiparticle density of 
states, and G has the value 3F&7/47r 3 ?i 3 . The equal- 
ity of the rates is evident upon substituting the explicit 
forms of the Fermi and Bose distribution functions. In 
the same manner, the emission and absorption rates are 
identical after the same substitutions and the change of 
variable E = E + v in P a . 

2.) The differences between the emission and absorp- 
tion rates and between the pair recombination and cre- 
ation rates have only a weak dependence on the lattice 
temperature as long as the electron temperature is near 
T c and the lattice temperature is low, say T p < T c /2. 
This conclusion is based on values for T c (7.75 meV) 
and A (19 meV) for YBCO. Wellstood, Urbina, and 
Clarke, 22 , assert that the difference between the emission 
rate and the absorption rate for a normal metal is equal 
to the difference between the rate electrons radiate to a 
zero temperature lattice and rate phonons radiate to a 
zero temperature electron gas. This result is only approx- 
imately valid in the gas of quasiparticles. The differences 
can be calculated from the previous pairs of equations 
by taking n(y) to be the phonon distribution function 
at T p . The difference between emission and absorption 
rates, for example, is: 



P e -P a = G dE 

Jo Jo 



E — A 



dvf e {E,E-u) 



■f/T P _ e u/T' 



The dependence on T p is contained in the last factor. 
For v of the order of A, T" of the order T c , and T in 
the range zero to T c /2, this factor only varies from 1.0 
to 0.915, showing therefore a weak dependence of the 
difference on T p . In an earlier calculation assuming s- 
wave symmetry of the order parameter, the differences 
were calculated explicitly for T = 0.8 T c and T p ranging 
from 0.1 T c to 0.5 T c . The difference varies less than 
10% for emission and absorption and less than 1% for the 
dominant creation and recombination. In view of these 
results, we ignore the dependence on T p and calculate r e 
on the basis of the radiation rate to a zero temperature 
lattice. A very similar argument applies to the difference 
between the pair recombination and pair creation rates. 

3.) Quasiparticle emission and absorbtion can only sat- 
isfy the energy and momentum conservation laws if the 
the quasiparticle velocity vpdE/de before emission or af- 
ter absorbtion is greater than the sound velocity This 
Cerenkov condition should be taken into account in the 
averaging near the Fermi surface that enters into calcula- 
tion of the electron-phonon spectral function a 2 F. In the 
integrals over quasiparticle energy above, the lower limit 
should be the energy E c at which quasiparticles reach 
the sound velocity rather than 0. The correction is of 
the order of the square of the ratio of the sound velocity 
to the Fermi velocity. Since s/vf << 1 for all supercon- 
ductors (s/vf ~ 1.5 x 10~ 2 in YBCO), the correction can 
be safely ignored. 



IV. PHONON LIFETIME EFFECTS 

In the above discussion, we have not distinguished the 
phonon temperature and the bath temperature. We now 
consider corrections arising from a more general treat- 
ment of the phonon distribution. The standard estimate 
of phonon mean free path for normal metals Hvp/kT 
gives 17 nm when the Fermi velocity is taken to be 2 x 10 5 
m/s, which is shorter than the 100 nm thickness of the 
experimental films. The estimate of Kaplan, Chu, Lan- 
genberg, Chang, Jafarey, and ScalapinoSi of quasipar- 
ticle and phonon lifetimes in an s-wave superconductor 
below the critical temperature gives a frequency- and 
temperature-dependent numerical factor of order unity 
multiplied by the characteristic time 



ph _ 



HN(a 2 



47T 2 iV(0) 



■A(0), 



(18) 



where N is the ion number density, (a 2 ) av is the average 
electron-phonon coupling constant, N(0) is the single- 
spin electronic density of states at the Fermi surface, and 
A(0) is the zero-temperature gap. Taking 2 ^ the values 
N = 13 per unit cell, (a 2 ) av = 5 meV, N(0) calculated 
from the free-electron theory with vf having the value 
quoted above, and A(0) = 19 meV, and converting the 
lifetime to a mean free path using the longitudinal sound 
velocity 4.2 x 10 3 m/s yields a path of the order of 10 3 
nm, an order of magnitude larger than the thickness of 
the experimental film. 

These estimates indicate that we are dealing with a 
case in which the phonon mean free path could be com- 
parable to the thickness of the sample. To deal with the 
general case, we follow Bezuglij and Shklovskij 13 , writing 
the kinetic equation for the phonon distribution function 
n(q, z) as 



dn{(j,z) 

dz 



n{cj,z)~n{T) 

Tph 



(19) 



where s z is the component of the sound velocity perpen- 
dicular to the plane of the film and n{T) is the ther- 
mal phonon distribution at the electron temperature. If 
phonons are reflected at the free surface of the film and 
transmitted with average coefficient a at the substrate in- 
terface, it is found that the phonon distribution function 
is a linear combination of two thermal distributions, one 
at the bath temperature and one at the electron temper- 
ature. The coefficients in the linear combination depend 
on the position within the film and on the direction of 
propagation of the phonons: 



n = A(z, B)n{T) + B(z, 6)n(T ). 



(20) 



The integrand in the expression for the energy transfer 
rate from quasiparticle gas to lattice contains the fac- 
tor n+1, while that for the reverse rate contains a fac- 
tor n. If Ea. (|20|l is substituted into these rates and ac- 
count is taken of the condition for equilibrium between 
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the lattice and the gas, the resulting rate contains a 
term with the factor 1 — A and a term with the factor 
Bn(To). For purposes of estimating the correction for 
finite phonon lifetime, we neglect the term proportional 
to n(To) compared to the 1 — A term on the ground that 
the phonon number is small at a temperature To which is 
much smaller than the Debye temperature. An estimate 
of the remaining term can be obtained by replacing 1 — A 
in the integral for the rate by its average value over the 
thickness of the film and over the directions of propaga- 
tion of the phonon. The remaining integral is the one we 
evaluated in the previous section. 

The resulting explicit expression for A is: 

A = , g / e-*/ l > , q z > 1 

1 - (1 - a)e- 2d ^ 1 e-( 2d -^ , q z < / - 

(21) 

Eq. |201 with [21 reflects the gradual change of the distri- 
bution from a nearly thermal distribution at the bath 
temperature at the substrate interface z — to an elec- 
tron temperature thermal distribution over the distance 
of a phonon mean free path. The transmission probabil- 
ity a can be determined in principle^ from the measured 
value of the thermal resistance of the film-substrate in- 
terface, defined as the ratio of AT at the interface to the 
product of the power dissipated per unit volume and the 
thickness of the film. The measured valu o 27 ' 28 for YBCO 
is about 1 x 1CP 3 Kcm 2 /W. The determination of a is 
affected by uncertainties due to the averaging and due 
to the temperature variation of the thermal resistance. 
A literal application of Eq.[14] of Ref.[2(| produces the 
average value 0.184 when d is of the order of or larger 
than I. When d <C I, Shklovskii shows that the effective 
a is 2d/l, which is 0.2 for the longest estimate of phonon 
mean free path above. We therefore accept 0.2 as a rea- 
sonable value. Sensitivity of the value of 1 — A to d/l and 
a are shown in Table [I] The longest estimate of phonon 
mean free path with the best estimate of the transmis- 
sion coefficient indicate that the energy transfer rate will 



be multiplied by a factor of 0.325 due to phonon lifetime 
effects. 



d/l 


a = 0.8 


a = 0.5 


a = 0.2 


6.0 


0.0333 


0.0208 


0.00833 


1.0 


0.190 


0.121 


0.0493 


0.1 


0.680 


0.553 


0.325 



TABLE I: Calculated values of the phonon lifetime factor 1— A 
for three values of the ratio of thickness d to phonon mean free 
path I and three values of the average transmission coefficient 
a. 



V. CONCLUSION 

At temperatures well below Tc, high electric fields and 
current densities can produce an instability that can be 
accounted for by a hot electron gas model in which the 
electronic temperature is elevated due to dissipation. The 
calculations presented here provide a quantitative justifi- 
cation for this scenario by showing that the temperature 
variation of the energy transfer rate between the lattice 
and the electrons is consistent with the position of the in- 
stability observed in YBCO films. They show also that if 
the phonon mean free path is not too small compared to 
the film thickness, the necessary temperature difference 
between electrons and lattice can be maintained. 
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